A systematic study of heavy quark energy loss in a hot and dense QCD medium 
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Within the framework of the Langevin equation, we perform a systematic study of heavy quark 
energy loss due to quasi-elastic multiple scatterings in a quark-gluon plasma created in relativistic 
heavy-ion collisions. We investigate how the initial configuration of the quark-gluon plasma as 
well as its properties affect the final state spectra and elliptic flow of D-mesons and non-photonic 
electrons. We find that both the geometric anisotropy of the initial quark-gluon plasma and the flow 
profiles of the hydrodynamic medium play important roles in the heavy quark energy loss process 
and the development of elliptic flow. The relative contributions of charm and bottom quarks is 
found to affect the transverse momentum dependence of the quenching and elliptic flow patterns 
of heavy flavor decay electrons; such influence depends on the interaction strength of heavy quarks 
with the medium. 



I. INTRODUCTION 

It has been well established that a deconfined state 
of QCD matter, called the quark-gluon plasma (QGP), 
can be created by colliding two heavy nuclei at ultra- 
relativistic energies, such as those at the Relativistic 
Heavy-Ion Collider (RHIC) and the Large Hadron Col- 
lider (LHC)[ll-Q. This highly excited form of matter 
shows properties similar to that of a nearly perfect fluid, 
such as strong collective flow which has been suc- 
cessfully described by relativistic hydrodynamic calcula- 
tions [li-El. 

Another interesting observation is the significant sup- 
pression of the production of high transverse momentum 
hadrons in high energy nucleus-nucleus collisions com- 
pared to that in binar y-sc aled proton-proton collisions at 
the same energies [17H19I ] . This suppression is commonly 
understood as the energy loss experienced by high en- 
ergy partons created in the initial hard collisions as they 
propagate through the deconfined hot and dense nuclear 
medium before fragmenting into hadrons [201 - I231 ] . Par- 
ton energy loss is thought to originate from a combina- 
tion of the collisional [25| and radiative 0, [2l| pro- 
cesses that jets experience when traversing the medium. 
Sophisticated jet quenching calculations have been per- 
formed for single inclusive hadron suppression [26l428l | as 
well as di-hadron }29l - l3l| and photon-hadron correlations 
f32Tj34j . Various energy loss schemes have been utilized in 
those phenomcnological applications, and for a detailed 
comparison between different formalisms, the reader is 
referred to Rcf.[35| and references therein. One of the 
goals of these calculations is the qualitative extraction of 
the jet transport parameters in the strongly-interacting 
QGP medium by comparing the calculations with the 
measured jet modification data. 

The production and the propagation of heavy quarks 
(charm and bottom) in relativistic heavy-ion collisions is 
of similarly great interest as they are expected to improve 
the constraints on our understanding of jet-medium in- 
teraction and provide further insights into the properties 
of the QGP. There has been observation of significant 



amount of suppression and elliptic flow for non-photonic 
electrons produced from the decay of heavy flavor mesons 
[36U38j ]. The suppression of heavy flavor decay elec- 
trons cannot be explained by the medium-induced ra- 
diation alone [39l442| . Heavy quarks, being massive, suf- 
fer medium-induced radiation suppressed in the forward 
cone compared to that for light partons; this is commonly 
referred to as the "dead cone effect" J43|. Thus, unlike 
the light partons whose energy loss is more dominated by 
medium-induced gluon radiation, the energy loss of heavy 
quarks is expected to be more dominated by collisional 
energy loss, especially at low to intermediate transverse 
momentum (4fl |43 - I47j . It has also been stated that 
the radiation in the backward region may have a larger 
impact on the heavy quark energy loss than expected 
before 0,111 . 

The large suppression of high momentum heavy quarks 
and their significant elliptic flow, both deduced from the 
measurement of non-photonic electrons, are usually re- 
garded as indication that heavy quarks may thermalize 
in the QGP jM H3, HIH ■ However, it was shown in 
Ref.[5Cjj that the thermalization time of charm quarks 
might be longer than the lifetime of the QGP phase for 
reasonable values of the heavy quark diffusion constant. 
Recent studies seem to indicate that the choice of a va- 
riety of medium-related parameters affects heavy quark 
spectra in relativistic heavy-ion collisions. For instance, 
Ref. 5l[ showed that with the identical transport coeffi- 
cients, different expansion scenarios can lead to a vari- 
ation by up to a factor of two for the suppression and 
elliptic flow of the heavy quark spectra. A similar inves- 
tigation of the medium flow effect on D-meson spectra 
was performed in Ref. (52j . In addition, uncertainties ex- 
ist in the calculation of initial heavy quark spectra, such 
as the charm-to-bottom quark ratio [39j , which may also 
affect the non-photonic electron spectra. 

In this paper, we carry out a systematic study of how 
heavy quark spectra in relativistic heavy-ion collisions 
depend on various ingredients in the phenomenological 
studies of heavy flavor energy loss, such as the initial 
production of heavy quarks, the geometry and the flow 
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properties of the hydrodynamic medium, and the cou- 
pling strength between heavy quarks and medium. We 
do not aim for a quantitative description of heavy flavor 
quenching data due to various uncertainties that we ex- 
plore here, such as the relative contributions from charm 
and bottom quarks to the production of non-photonic 
electrons. The measurements of the suppression and 
flow for charm and bottom mesons separately will greatly 
improve this situation and put more constraints on the 
interpretation of the data. For our purpose, we shall 
only consider the collisional energy loss experienced by 
heavy quarks propagating through the QGP medium. 
Our study should be applicable to heavy quark energy 
loss in the region from low to intermediate transverse 
momenta. The contribution from radiative energy loss 
to heavy quark evolution in medium will be included in 
an upcoming effort. 

In the limit of multiple quasi-elastic interactions where 
the energy transfer during each collision is small, the 
propagation of heavy quarks through a QGP medium can 
be treated as Brownian motion. We will follow Ref.[50l| 
and study heavy quark evolution in the medium in the 
framework of Langevin equation [3, l5ll - [5a |. In this 
framework, one of the crucial parameters is the diffu- 
sion coefficient, which has been widely discussed from a 
microscopic point of view (57l - [60| . 

The paper is organized as follows. In Sec. II, we 
present the methodology utilized for simulating heavy 
quark evolution in a dynamic medium as produced in 
rclativistic heavy-ion collisions. In Sec. Ill, we present 
the numerical results of the nuclear modification factor 
and elliptic flow for heavy quarks, heavy flavor mesons, 
and non-photonic electrons, where their dependence on 
various inputs are discussed in details. A summary and 
outlook is presented in Sec. IV. 



II. METHODOLOGY 

In this work, we only consider the collisional energy 
loss experienced by heavy quarks propagating in the 
QGP medium. The multiple quasi-elastic scatterings of 
heavy quarks inside a thermalized medium can be treated 
as brownian motion and be described by the Langevin 
equation: 



(1) 



In principle, the noise term £ may depend on the momen- 
tum of the particle, but we assume for simplicity that this 
is not the case here. The random momentum kicks satisfy 
the following correlation relation: 



(2) 



adopted, 

p(t + At) = p(t) - d Ito (p(t))At + |*At, 
^(t)^(t-nAt)) = ^S 0n , (3) 

where the force due to the drag is given by 

dito(p) = Vd(p)p- (4) 

Assuming the energy transfer is small, it can be shown 
that the fluctuation-dissipation relation applies, which 
indicates: 



Vd(p) 



2TE 



(5) 



Furthermore, according to Eq. (J3j , Gaussian noise with 
width r = \J nj At will be used to generate the random 
momentum kicks in our calculation. The diffusion coef- 
ficient is related to the drag term via: 
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In the simulation of the momentum evolution of heavy 
quarks from Langevin equation, the Ito discretization is 



The detailed information about the produced highly 
excited medium in rclativistic heavy-ion collisions is es- 
sential to heavy quark energy loss and the development 
of elliptic flow. Throughout our calculation, the QGP 
medium is described by the fully (3+l)-D relativistic 
ideal hydrodynamic model developed in Ref. [l2| . In our 
work, two different initial condition models, a Glauber 
IH IIS] as well as a KLN-CGC [i| model, are utilized 
for the description of the initial energy distribution of the 
medium before the hydrodynamic evolution. The com- 
parison between the two will allow for a study of the 
sensitivity of heavy-quark observables to the initial spa- 
tial make-up of the system. We will focus on mid-central 
Au-Au collisions at RHIC with a center-of-mass energy 
Y / J/vw=200 GeV per nucleon pair and take the impact 
parameter b = 6.5 fm throughout the calculation. 

One of the key assumptions of the hydrodynamic cal- 
culation is that the QCD medium experiences a sudden 
thermalization to form a QGP matter at an initial time 
ro (chosen to be 0.6 fm/c in this work) at which the hy- 
drodynamic evolution commences. The pre-equilibrium 
phase of the medium cannot be modeled via hydrody- 
namics - here we assume that heavy quarks stream freely 
prior to QGP formation. Any heavy quark leaving the 
QGP will stream freely as well. 

The hydrodynamic simulation provides us with the 
time-evolution and the spatial distribution of temper- 
ature and flow velocities for the thermalized QGP 
medium. In such a dynamic medium, heavy quark evo- 
lution is treated as follows: for every Langevin time step 
we boost the heavy quark to the local rest frame of the 
fluid cell through which it propagates. The Langevin ap- 
proach is then applied to obtain the momentum evolution 
of heavy quark. After that we boost back to the global 
computation frame. 
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Since the initial production of heavy quarks is domi- 
nated by the processes with high transverse momentum 
transfer, perturbative QCD is applied to model the ini- 
tial momentum distribution of the heavy quarks, prior to 
their propagation through the QGP medium. For sim- 
plicity, we sample the initial transverse momentum of 
heavy quarks according to the following power-law dis- 
tribution [in, 



1 



dN 



d 2 p T (p| + A2 ) Q ' 



(7) 



By comparing this expression to leading-order perturba- 
tive QCD calculations, we fix the above parameters a 
and A to be: a = 3.9 and A = 2.1 for charm quarks, 
and a = 4.9 and A = 7.5 for bottom quarks. In this 
work, we focus on the energy loss of heavy quarks at 
mid-rapidity and therefore assign no initial longitudinal 
momentum to the heavy quarks in the initial state. We 
have checked that the introduction of initial longitudinal 
momenta that are uniform around the mid-rapidity re- 
gion (—1 < rj < 1) does not affect our final transverse 
momentum spectra and does not affect the systematics 
we are about to explore. The relative normalization (ra- 
tio) of charm and bottom quarks is not fixed, but rather 
serves as a free parameter in our simulation. Later we will 
investigate the effect of this normalization on the quench- 
ing and the elliptic flow of heavy flavor decay electrons. 

The initial spatial distribution of heavy quarks in the 
transverse plane is sampled according to the distribu- 
tion of binary collisions as calculated from a Monte- Carlo 
Glauber model. With the spatial and momentum ini- 
tialization of heavy quarks, we are able to simulate their 
time evolution inside the QGP medium in the framework 
of Langevin equation as described above. After passing 
through the medium, their fragmentation into heavy fla- 
vor mesons and the subsequent decay into electrons are 
simulated via Pythia 6.4 [64j. In the end, the final state 
particles in the mid-rapidity region (— 1 < rj < 1) are 
selected and their momentum distribution are utilized to 
calculate the elliptic flow coefficient «2 and the nuclear 
modification factor Raa as follows: 



v 2 (p T ) = (cos(2ct>)) = 

\ Px 



Raa 



(dN/d PT ) Rn 
{dN/dp T \ 



(8) 



init 



Note that when heavy quarks are directly analyzed, the 
denominator and the numerator of Raa are the initial 
heavy quark distribution and the distribution of those 
surviving from the energy loss and passing through the 
medium. When analyzing heavy flavor mesons or elec- 
trons, the denominator represents the spectra of the cor- 
responding particles fragmented/decayed directly from 
the initial heavy quarks, while the numerator represents 
those produced from the heavy quarks after transporting 
through the QGP medium. 



III. NUMERICAL RESULTS 

In this section, we will present the numerical results 
for the simulation on heavy quark evolution in the QGP 
medium produced by Au+Au collisions at RHIC energy. 
We will investigate the impact various parameter choices 
of the calculation have on the final spectra and elliptic 
flow of the heavy quarks, heavy mesons and their decay 
electrons. 



A. Charm Quark Energy Loss and Flow 

The energy loss of heavy quarks and the develop- 
ment of elliptic flow crucially depend on the geometri- 
cal shape and dynamical evolution of the thermalized 
QGP medium that the heavy quarks traverse. The total 
energy loss of heavy quarks is mostly controlled by the 
overall magnitude of the energy density of the medium, 
while the elliptic flow is more sensitive to the geome- 
try of the medium as it characterizes the anisotropy of 
the final transverse momentum spectra. In typical non- 
central nucleus-nucleus collisions, the overlap region of 
the two nuclei is anisotropic in the transverse plane, thus 
resulting in the anisotropy of the produced hot and dense 
medium. Due to the different pressure gradients in differ- 
ent directions, different radial flows are built up during 
the hydrodynamic evolution of the thermalized QGP. 

In such an anisotropic dynamical medium, there ex- 
ist two factors affecting the anisotropy of heavy quark 
energy loss: the different path lengths and the differ- 
ent flow profiles experienced by the heavy quarks travel- 
ing in different directions. Longer paths will be traveled 
through by heavy quarks moving in the out-of-plane (y) 
direction than in the in-plane (x) direction, where the 
reaction plane is defined to be spanned by the impact 
parameter and the beam axis directions. Thus in ab- 
sence of collective flow for the medium, heavy quarks, 
after passing through such anisotropic medium, would 
have larger momentum in the x direction than in the y 
direction, > {p y }, resulting in a positive elliptic flow. 
In addition, the collective flow of the medium also con- 
tributes positively to heavy quark elliptic flow since the 
push of the radial flow is more prominent in the x direc- 
tion. Therefore, the total elliptic flow developed during 
the propagation of heavy quarks in such an anisotropic 
hydrodynamic medium is due to a combination of these 
two factors. 

We can separate these two effects in the simulation by 
turning on or off the coupling of the collective flow of the 
thermalized medium to the evolving heavy quarks. The 
decoupling from the collective flow can be accomplished 
by not boosting the heavy quarks into the respective rest 
frame of the fluid cell for the Langevin evolution. The 
comparison between the heavy quark evolution with and 
without coupling to the collective flow is shown in FigfTJ 
where the upper plot shows the nuclear modification fac- 
tor Raa and the lower plot shows the elliptic flow vi of 
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FIG. 1: (Color online) A comparison between the influence 
of QGP media with and without collective flow on (a) Raa 
and (b) V2 of charm quarks. Both media are generated with 
Glauber initial condition. 



the charm quarks as a function of the transverse momen- 
tum. We show results for two different values of diffusion 
coefficient D = 1.5/(2ttT) and D = 6/(2ttT). 

The effect of the collective flow of the medium on the 
heavy quark energy loss can be clearly seen from the 
plot of the nuclear modification factor Raa (FigHJ . It is 
negligible at high transverse momenta, and becomes ob- 
servable at intermediate transverse momentum regime. 
Due to the push by the radial flow, heavy quarks are 
less suppressed (i.e. have a larger Raa) at larger trans- 
verse momenta, since the radial flow effectively trans- 
ports low momentum heavy quarks to larger transverse 
momenta. Similar effects stemming from the elliptic flow 
of the medium are observed for the heavy quark ellip- 
tic flow coefficient V2- At low transverse momenta, the 



collective flow of the medium presents a significant in- 
fluence on the charm i>2. At high transverse momenta, 
the collective flow effect is small, thus the development 
of charm V2 is dominated by the geometric anisotropy of 
the medium. The dominance of the medium collective 
flow at low transverse momentum for vi might indicate 
that low transverse momentum charm quarks are more 
likely to lose a significant amount of their momenta and 
thcrmalizc in the medium, and thus flow more like the 
thermalizcd medium. 

A closer observation suggests that with a decrease 
of the diffusion coefficient, i.e., an increase of the cou- 
pling strength, the influence of the geometric asymme- 
try becomes more dominant. For instance, Fig[T] reveals 
that for D = 6/(2irT), the geometric asymmetry of the 
medium contributes to only approximately half of charm 
i>2 at the peak value (around pr = 1-5 GeV). How- 
ever, for D = 1.5/(27rT), such contribution increases to 
more than 80% at the corresponding peak value (around 
Pt = 3 GeV) . Note that such an increase of the geometric 
contribution is not unlimited. With further reduction of 
the value of the diffusion coefficient (D < 1.5/(27rT)), 
i.e., a larger coupling between heavy quarks and the 
medium, the energy loss of charm quarks will be so in- 
tense that all of them will be captured by the medium. 
In that limit, charm quarks thermalizc with the medium 
during the QGP lifetime [50| . and therefore, their v-i will 
entirely follow collective flow of the medium. In our sim- 
ulation, the choice of D ~ 1.5/ (2ttT) provides the largest 
elliptic flow for the final heavy quarks. 

We may further investigate the effect of the spatial 
medium distribution on the heavy quark energy loss and 
the development of heavy quark elliptic flow by utilizing 
different initial conditions for the hydrodynamic simu- 
lation of the QGP. Two different initial condition mod- 
els are widely used for the initialization of the energy 
density distribution prior to the hydrodynamic evolution: 
the Glauber model [6l|, |62[ and KLN parametrization of 
the Color Glass Condensate (CGC) model (63[. These 
two models provide initial energy density profiles with 
different anisotropics in the transverse plane. In partic- 
ular, the KLN-CGC model exhibits a larger eccentricity 
£2 = (y 2 — x 2 ) / (y 2 + x 2 ) than the Glauber model, which 
will manifest itself in larger elliptic flow coefficients for 
the heavy quarks. 

The comparison between these two initial condition 
models is shown in Fig|2] where the upper frame of the 
figure shows the nuclear modification factor Raa and 
lower frame shows the elliptic flow V2- As expected, a 
significantly larger elliptic flow is observed for the charm 
quarks traveling through the hydrodynamic medium with 
KLN-CGC initial condition than those with Glauber ini- 
tial condition. As indicated by Fig|3J the difference 
can be as large as 20% for D = 6/(2ttT) and 40% for 
D = 1.5/(27rT). We also observe that while v-i is sen- 
sitive to the choice of the initial condition, the nuclear 
modification factor Raa is not significantly affected by 
the choice of these two hydro initial conditions. This is 
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FIG. 2: (Color online) A comparison between the influence 
of QGP media with the Glauber and the KLN-CGC initial 
conditions on (a) Raa and (b) vi of charm quarks. 



FIG. 3: (Color online) (a) Raa and (b) vi of Do mesons. 
The QGP medium is generated with the KLN-CGC initial 
condition. 



due to Raa being controlled by the overall normaliza- 
tion of the density profile in the hydrodynamic evolution 
which has been tuned to describe the properties of bulk 
matter, such as the it and K spectra. 



B. D Mesons and Heavy Decay Electrons 

In the above discussion, we have focused on the effects 
of initial conditions and medium parameters on heavy 
quark energy loss and the development of heavy quark 
elliptic flow. Now we investigate the corresponding sensi- 
tivities for heavy flavor mesons and their decay electrons. 
Since the KLN-CGC initial condition provides a larger 
eccentricity for the initial energy density profile and thus 



produces a larger elliptic flow in the heavy quarks during 
their medium evolution, we shall use it for the remain- 
der of our analysis. This is merely to obtain the largest 
possible values for the final elliptic flow, since most of 
the previous calculations seem to under-predict the ellip- 
tic flow data for non-photonic electrons once the model 
parameters have been tuned to describe the measured 
nuclear modification factor. 

FigJ3]and Fig|5]display the numerical results of the nu- 
clear modification factor Raa (upper panel) and elliptic 
flow V2 (lower panel) for D mesons and £>-decay elec- 
trons. Three different values of diffusion coefficients are 
used for comparison D = 1.5/(2ttT), D = 3/(2ttT), and 
D = 6/(2%T). We observe that the transverse momen- 
tum dependence of Raa and i>2 are similar to that for 
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FIG. 4: (Color online) (a) Raa and (b) V2 of electrons decayed FIG. 5: (Color online) A comparison of (a) Raa and (b) V2 
from charm quarks. The QGP medium is generated with the between charm quarks, Do mesons and electrons. The QGP 
KLN-CGC initial condition. medium is generated with the KLN-CGC initial condition. 



charm quarks as shown in the previous figures. 

Fig (5] provides a more direct comparison of the nuclear 
modification of Raa (upper frame) and elliptic flow 
(lower frame) between charm quarks, D meson and D- 
decay electrons. Here the result is shown for a diffusion 
coefficient of D = 1.5/ (2irT); other values of the diffu- 
sion coefficient display similar systcmatics. One observes 
that all the nuclear modification factors Raa curves de- 
crease with increasing transverse momentum, and satu- 
rate above pt ~ 3 — 4 GeV, while the elliptic flow curves 
first increase, reach a peak and then decrease (and sat- 
urate). The peak values of V2 is the largest for charm 
quarks and the smallest for electrons. This decrease of 
the momentum anisotropy can be understood as a result 
of additional randomization of the momentum space dur- 
ing the hadronization and decay processes. Moreover, the 



values of the transverse momentum at which V2 reaches 
the peak value shift to lower regimes from charm quarks 
to D mesons and to D-decay electrons. This pattern is 
caused by the decrease of the momentum from the par- 
ent particles to the daughter particles during heavy quark 
fragmentation and heavy flavor meson decay processes. 

For the heavy flavor decay electron spectra, another 
important factor is the relative contributions from charm 
vs. bottom quarks. Since charm and bottom quarks 
have different masses, they are produced with different 
initial transverse momentum distributions, and experi- 
ence different energy loss and coupling to the collective 
flow in medium. This manifests itself in different Raa 
and vi systematics for D and B mesons respectively and 
subsequently translates into different behavior for their 
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FIG. 6: (Color online) A comparison of (a) Raa and 
(b) V2 of non-photonic electrons between different initial 
charm/bottom ratios. Set D = 6/(2-7rT) and the QGP 
medium is generated with the KLN-CGC initial condition. 



FIG. 7: (Color online) A comparison of (a) Raa and 
(b) V2 of non-photonic electrons between different initial 
charm/bottom ratios. Set D = 1.5/(27rT) and the QGP 
medium is generated with the KLN-CGC initial condition. 



respective decay electrons. The electrons at lower trans- 
verse momentum arc dominated by charm quark (i.e. D 
meson) decay, while in the high transverse momentum 
regime bottom quark dominate as sources for these de- 
cay electrons. Since there are multiple uncertainties af- 
fecting the relative normalization of charm and bottom 
quark production, for example the scale dependence in 
the perturbative QCD calculation of initial heavy quark 
production [39j . we treat the ratio of charm and bottom 
quarks as a free parameter for our calculation, and in- 
vestigate how the variation of this ratio affects the final 
non-photonic electron distributions. 

The results are shown in FigJB] and FigiT] for two dif- 
ferent values of diffusion coefficients, D = 1.5/(2irT) 
and D = 6/(2ttT), respectively. In each figure, the up- 



per frame shows the nuclear modification factor Raa 
and the lower frame shows the elliptic flow «2- Wc 
compare four different initializations here - pure charm, 
pure bottom, and two mixtures of charm and bottom 
quarks: 99.2% charm quarks with 0.8% bottom quarks, 
and 98.5% charm quarks with 1.5% bottom quarks. As 
shown in [39| , the bottom quark contribution to the elec- 
tron spectra may start dominating over charm quark con- 
tribution at transverse momentum as low as 3 GeV or as 
high as 9 GeV. Our two hybrid mixtures of charm and 
bottom quarks have about a factor of 2 difference in their 
ratio, representing an estimate of the uncertainties due 
to our limited control of the proton-proton baseline. 

One observes from these two figures that the nuclear 
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modification factor Raa and the elliptic flow of heavy 
flavor decay electrons are very different for the pure 
charm vs. pure bottom scenario. Bottom quarks are less 
suppressed than charm quarks at high transverse mo- 
menta, thus less enhancement is observed at low trans- 
verse momenta in the Raa- The magnitude of the el- 
liptic flow coefficient %)% is much smaller for decay elec- 
trons stemming from bottom decay than from charm de- 
cay, again due to the reduced energy loss experienced 
by the bottom quarks. However, we also observe a dif- 
ference in the transverse momentum dependence: while 
the elliptic flow coefficient v-x of electrons from charm 
decays has a peak value at intermediate transverse mo- 
mentum, that for bottom decays increases monotonically 
with increasing transverse momentum (and then satu- 
rates). These result from the fact that charm quarks 
contribute mostly to the production of low transverse 
momentum non-photonic electrons while bottom quarks 
contribute mostly to high transverse momentum elec- 
trons. 

Due to the different behavior of charm vs. bottom 
decay electrons, the electrons from a mixture of charm 
and bottom decays exhibit a very rich structure. Both 
Raa and vi trend similar to the pure charm initializa- 
tion at low transverse momentum and converge to the 
values of the pure bottom quark scenario at high trans- 
verse momenta. In the intermediate transverse momen- 
tum region where the transition from charm dominance 
to bottom dominance in the origin of the decay electrons 
takes palce, a non-monotonic transverse momentum de- 
pendence of Raa and vi is observed: a dip-peak structure 
for Raa and a peak-dip structure for r>2- Such a non- 
monotonic behavior is more prominent for the smaller 
value of the diffusion coefficient D = 1.5/(2ttT) (FigJTJ), 
since a smaller value of the diffusion coefficient increases 
the interaction with the medium and thus the energy loss 
of charm quarks and their elliptic flow, while such an en- 
hancement is far less for bottom quarks due to their larger 
mass. 

Current experimental results seem not able to deter- 
mine whether such a peak-dip structure is present or not 
in the non-photonic electron elliptic flow v-i due to large 
experimental error bars. Further improvement of the 
measurement of the detailed transverse momentum de- 
pendence of non-photonic electrons would be helpful for 
the determination of the diffusion coefficient and there- 
fore the coupling strength between heavy quarks and the 
QGP medium. 

Another important effect seen in Fig[B] and FigJT] is 
the significant sensitivity of heavy flavor decay electron 
t>2 to the initial charm-to-bottom quark ratio. For in- 
stance, a 0.7% difference in the mixing ratio between 
charm and bottom quarks in our simulation leads to a 
variation of approximately 25% in for a diffusion coef- 
ficient of D = 6/(2nT) and over 30% for D = 1.5/(2ttT). 
As has been discussed earlier, significant uncertainties re- 
garding the initial heavy quark spectra are still present 
in our current phenomcnological calculations, and thus 



provide a sizable uncertainty for the prediction of the 
quenching and elliptic flow of non-photonic electrons. 



IV. SUMMARY AND OUTLOOK 

In this work, we have performed a systematic study 
of heavy quark energy loss in a hot and dense QGP 
medium produced in rclativistic heavy-ion collisions. The 
Langcvin approach is utilized to simulate the heavy quark 
evolution in the medium due to quasi-elastic multiple 
scatterings. Numerical results are presented for both 
the nuclear modification factor and the elliptic flow of 
heavy quarks, heavy flavor mesons and their correspond- 
ing non-photonic decay electrons. We have investigated 
in details how the final nuclear modification factor and 
elliptic flow are affected by various components of the 
model, such as the geometry and the collective flow of 
the hydrodynamic medium, the initial production ratio 
of charm to bottom quarks and the coupling strength 
between the heavy quarks and the medium. 

We have focused on two particular properties of the 
medium that affect the heavy quark energy loss - its ge- 
ometric anisotropy and its collective flow. It is found that 
the geometric anisotropy dominates the final heavy quark 
distributions in the high transverse momentum region, 
while the collective flow of the medium dominates the 
low momentum region. The impact of the QGP geome- 
try on the heavy quark energy loss has been explored by 
comparing Glauber and KLN-CGC initializations for the 
hydrodynamic medium. We found that while a similar 
nuclear modification factor Raa is observed for both ini- 
tial condition models, a significantly higher heavy quark 
elliptic flow V2 is found for the KLN-CGC initial condi- 
tions. We have further investigated the sensitivity of the 
spectra and elliptic flow of non-photonic electrons to the 
relative contributions from charm and bottom quarks. 
It is found that a less than 1% difference in the initial 
charm-to-bottom ratio can lead to more than 30% vari- 
ation of the non-photonic electron spectra. Therefore, 
narrowing down these uncertainties is essential for a bet- 
ter understanding of the interaction dynamics between 
heavy quarks and the QGP medium. 

Though we have not preformed a detailed quantitative 
comparison to the experimental measurements, most of 
our results can be tested with future measurements. For 
instance, the spectra of the heavy flavor decayed elec- 
trons show a dip-peak structure of Raa and a peak-dip 
structure of V2 (see Fig[7]) for large values of the quark- 
medium coupling; such structures would disappear if the 
coupling becomes weaker (Fig (5]). Due to the large error 
bars of the data in the high transverse momentum region, 
current experimental observations can not yet distinguish 
between these two scenarios. Future high precision mea- 
surements that can be compared to this prediction will 
help determining the coupling strength between heavy 
quarks and the QGP medium. 

While our study constitutes an important move for- 
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ward to quantitatively understand the interaction dy- 
namics of heavy quarks in a hot and dense QGP, it can 
be further improved in several directions, which we leave 
for future work. Here, we have utilized the Langevin 
approach for the simulation of heavy quark evolution 
as affected by multiple quasi-elastic scatterings in the 
QGP medium, in which many details about the micro- 
scopic structures of the interaction between heavy quarks 
and the medium are encoded in the diffusion coefficient. 
However, the heavy quark scattering cross section in the 
non-perturbative regime remains largely unknown, re- 
sulting in different predictions for this transport coeffi- 
cient. In addition, most current calculations do not ac- 
count for any interaction in the pre-equilibrium phase of 
the medium prior to QGP formation. However, turbu- 
lent chromo-electromagnetic fields in this pre-equilibrium 
stage may have an influence on the heavy quark motion 
that cannot be overlooked (65|. Another interesting ef- 
fect on the heavy quark evolution in medium is the ra- 



diative energy loss induced by multiple scatterings. This 
effect has been discussed in other frameworks like the 
Boltzmann equation [gf|, but is still absent in the cur- 
rent Langevin algorithm and may be important, espe- 
cially in the high transverse momentum regime. In a 
follow-up work, we shall address some of these questions, 
in particular how the pre-equilibrium phase of the QGP 
evolution and medium-induced gluon radiation affect the 
heavy quark energy loss. 
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